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Abstract 

Production of exclusive dijets in diffractive deep inelastic e^p scattering has been 
measured with the ZEUS detector at HERA using an integrated luminosity of 
372 pb“^. The measurement was performed for 7 *-p centre-of-mass energies in the 
range 90 < VU < 250 GeV and for photon virtualities > 25 GeV^. Energy flows 
around the jet axis are presented. The cross section is presented as a function of (5 
and (/>, where /3 = x/xp, x is the Bjorken variable and xp is the proton fractional 
longitudinal momentum loss. The angle 4> is defined by the 7 *-dijet plane and the 
7 *-e^ plane in the rest frame of the diffractive final state. The (j) cross section is 
measured in bins of (5. The results are compared to predictions from models based 
on different assumptions about the nature of the diffractive exchange. 
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1 Introduction 


The first evidence for exclusive dijet production at high-energy hadron colliders was 
provided by the CDF experiment at the Fermilab Tevatron pp collider [T] and had an 
important impact on theoretical calculations of exclusive Higgs boson production at the 
Large Hadron Collider. This p^er describes the first measurement of exclusive dijet pro¬ 
duction in high energy electroitJ-proton scattering. A quantitative understanding of the 
production of exclusive dijets in lepton-hadron scattering can improve the understanding 
of more complicated processes like the exclusive production of dijets in hadron-hadron 
scattering [2] or in lepton-ion scattering at a future eRHIC accelerator |3]. 

A schematic view of the diffractive production of exclusive dijets, e+p —)■ e-l-jetl-|-jet2-|-p, 
is shown in Fig. [T] In this picture, electron-proton deep inelastic scattering (DIS) is 
described in terms of an interaction between the virtual photon, 7*, and the proton, which 
is mediated by the exchange of a colourless object called the Pomeron (P). This process 
in the 7*-P centre-of-mass frame is presented in Fig. HI where the lepton and jet planes 
are marked. The lepton plane is defined by the incoming and scattered electron momenta. 
The jet plane is defined by the jet momenta, which are always back-to-back, and the 
virtual photon momentum. The angle between these planes is labelled 0. The jet polar 
angle is defined with respect to the virtual photon momentum and called 6. 

The production of exclusive dijets in DIS is sensitive to the nature of the object exchanged 
between the virtual photon and the proton. Calculations of the single-differential cross 
section of dijet production as a function of 0 in fct-factorisation |1] and collinear factor¬ 
isation in] have shown that, when the quark and antiquark jets are indistinguishable, the 
cross section is proportional to 1 -|- A(pr,jet) cos 20, where PTjet is the jet transverse mo¬ 
mentum. It was pointed out for the first time by Bartels et al. mu that the parameter 
A is positive if the quark-antiquark pair is produced via the interaction of a single gluon 
with the virtual photon and negative if a system of two gluons takes part in the inter¬ 
action. The absolute value of the A parameter is expected to increase as the transverse 
momentum of the jet increases. 

The production of exclusive dijets is also sensitive to the gluon distribution in the proton 
and is a promising reaction to probe the off-diagonal (generalised [7|) gluon distribution. 
The off-diagonal calculations predict a larger cross section compared to calculations based 
on conventional gluon distributions. In this context, the exclusive production of dijets is 
a complementary process to the exclusive production of vector mesons which has been 
extensively studied at HERA [Hl-fTb]. 

This paper describes the measurement of differential cross sections as a function of f3 and 
in bins of /9 as a function of 0. The former quantity is defined as 0 = x/xp, where x 
is the Bjorken variable and xp is the fractional loss of proton longitudinal momentum. 

^ Here and in the following the term “electron” denotes generically both the electron and the positron. 
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The results of this analysis are compared to predictions from the Two-Gluon-Exchange 
model and the Resolved-Pomeron model of Ingelman and Schlein |18] . 


2 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [12]. A brief outline 
of the components that are most relevant for this analysis is given below. 

In the kinematic range of the analysis, charged particles were tracked in the central tracking 
detector (CTD) |20| and the microvertex detector (MVD) [2T| . These components operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angl^ region 15° < 6 < 164°. The MVD silicon tracker consisted of a barrel (BMVD) 
and a forward (FMVD) section. The BMVD contained three layers and provided polar- 
angle coverage for tracks from 30° to 150°. The four-layer FMVD extended the polar-angle 
coverage in the forward region to 7°. After alignment, the single-hit resolution of the MVD 
was 24 pm. The transverse distance of closest approach (DCA) of tracks to the nominal 
vertex in X-Y was measured to have a resolution, averaged over the azimuthal angle, of 
(46 © 122/pr) pni, with px in GeV. For CTD-MVD tracks that pass through all nine CTD 
superlayers, the momentum resolution was a{pT)/pT = 0.0029pT © 0.0081 © 0.0012 /pt, 
with pt in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [22| consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
was subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a{E)/E = 0.18/\/E for electrons and 
a{E)/E = 0.35/\/E for hadrons, with E in GeV. 

The position of electrons scattered at small angles to the electron beam direction was 
determined with the help of RHES [23], which consisted of a layer of approximately 10 000 
(2.96 X 3.32 cm^) silicon-pad detectors inserted in the RCAL at a depth of 3.3 radiation 
lengths. 

The luminosity was measured using the Bethe-Heitler reaction ep ejp by a luminosity 
detector which consisted of independent lead-scintillator calorimeter [23] and magnetic 
spectrometer [23] systems. The fractional systematic uncertainty on the measured lumin¬ 
osity was 2% [26] . 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
nominal proton beam direction, referred to as the “forward direction”, and the X axis pointing towards 
the centre of HERA. The coordinate origin is at the centre of the CTD. The pseudorapidity is defined 
as ?7 = — In (tan |), where the polar angle, 0, is measured with respect to the Z axis. 
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3 Monte Carlo simulation 


Samples of Monte Carlo (MC) events were generated to determine the response of the 
detector to jets of hadrons and the correction factors necessary to obtain the hadron-level 
jet cross sections. The hadron level is defined in terms of hadrons with lifetime > 10 ps. 
The generated events were passed through the GEANT 3.21-based [27] ZEUS detector- 
and trigger-simulation programs [19]. They were reconstructed and analysed by the same 
program chain used for real data. 

In this analysis, the model SATRAP [281129] as implemented in the RAPGAP [30] program 
was used to generate diffractive events. SATRAP is a colour-dipole model [31] which 
includes saturation effects. It describes DIS as a fluctuation of the virtual photon into a 
quark-antiquark dipole which scatters off the proton. The GTEQ5D [32] parameterisation 
was used to describe the proton structure. Hadronisation was simulated with the JETSET 
7.4 [331134] program which is based on the Lund string model [33]. Radiative corrections 
for initial- and final-state electromagnetic radiation were taken into account with the 
HERAGLES 4.6.6 [361137] program. The diffractive MG was weighted in order to describe 
the measured distributions (see Section [3]) • 

The proton-dissociation process was modelled using the EPSOFT [381139] generator. The 
production of dijets is not implemented in EPSOFT. Therefore dijets with proton disso¬ 
ciation were simulated with SATRAP, where the intact proton was replaced with a dis¬ 
sociated proton. Such a solution is based on the factorisation hypothesis which assumes 
that the interaction at the lepton and at the proton vertex factorises. The factorisation 
hypothesis has been verified for diffractive processes in ep collisions at HERA [401145] . 

To estimate the non-diffractive DIS background, a sample of events was generated using 
HERAGLES 4.6.6 [301137] with DJANGOH 1.6 [44] interfaced to the hadronisation pro¬ 
cess. The QGD cascade was simulated using the colour-dipole model (GDM) [451147] as 
implemented in ARIADNE 4.08 [451140] . 

To estimate the background of diffractive dijet photoproduction, a sample of events was 
generated using the PYTHIA 6.2 [50] program with the GTEQ4L [51] parton density 
function of the proton. The hadronisation process was simulated with JETSET 7.4. 

For the model predictions, events were generated using RAPGAP where both the Resolved- 
Pomeron model and the Two-Gluon-Exchange model are implemented. The hadronisation 
was simulated with ARIADNE. The generated events do not include proton dissociation. 

In this analysis, the number of diffractive MG events was normalised to the number of 
events observed in the data after all selection cuts and after subtraction of background 
from photoproduction and non-diffractive DIS. The numbers of background events were 
estimated based on generator cross sections. 
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4 Event selection and reconstruction 


This analysis is based on data collected with the ZEUS detector at the HERA collider 
dnring the 2003-2007 data-taking period, when electrons or positrons of 27.5 GeV were 
collided with protons of 920 GeV at a centre-of-niass energy of ^/s = 318 GeV. The data 
sample corresponds to an integrated luminosity of 372 pb“^. 

A three-level trigger system was used to select events online |19ll52] . At the first level, 
only coarse calorimeter and tracking information were available. Events consistent with 
diffractive DIS were selected using criteria based on the energy and transverse energy 
measured in the GAL. At the second level, charged-particle tracks were reconstructed 
online by the ZEUS global tracking trigger [531 l5l] . which combined information from the 
GTD and MVD. These online tracks were used to reconstruct the interaction vertex and 
to reject non-ep background. At the third level, neutral current DIS events were accepted 
on the basis of the identification of a scattered electron candidate using localised energy 
depositions in the GAL. 

The scattered electron was identified using a neural-network algorithm [SS]- The recon¬ 
struction of the scattered electron variables was based on the information from the GAL. 
The energy of electrons hitting the RGAL was corrected for the presence of dead material 
using the rear presampler detector [SB]. Energy-flow objects (EFOs [FTIISH] ) were used to 
combine the information from the GAL and the GTD. 

4.1 DIS selection 

A clean sample of DIS events with a well-reconstructed electron was selected by the fol¬ 
lowing criteria: 

• the electron candidate was reconstructed with calorimeter information and was re¬ 
quired to have energy reconstructed with double-angle method [59], U' > 10 GeV 
and, if reconstructed in the GTD acceptance region, also an associated track; 

• the reconstructed position of the electron candidate in the GAL was required to be 
outside the regions of GAL in which the scattered electron might have crossed a 
substantial amount of inactive material or regions with poor acceptance; 

• the vertex position along the beam axis was required to be in the range \ Z^t^\ < 30 cm; 

• T^had/A'tot > 0.06, where Uhad is the energy deposited in the hadronic part of the GAL 
and E^ot is the total energy in the GAL; this cut removes purely electromagnetic 
events; 

• 45 < (U — Pz) < 70 GeV, where E is the total energy, E = ^. Ei, Pz = and 

pz^i = Ei cos 6i, where the sums run over all EFOs including the electron; this cut 
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removes events with large initial-state radiation and further reduces the background 
from photoproduction. 

Events were accepted if > 25 GeV^ and 90 < hh < 250 GeV. In this analysis, the 
photon virtuality, Q^, and the total energy in the virtual-photon-proton system, W, were 
reconstructed using the double-angle method which was found to be more precise than 
other reconstruction methods in the kinematic region of this measurement j60]. The 
inelasticity, y, which was reconstructed with the electron method was limited to the range 
0.1 < y < 0.64. The limits come from the selection criteria applied to other variables 
reconstructed with the double-angle method. The use of two methods to reconstruct DIS 
kinematic quantities, increases the purity if the sample. 

4.2 Diffractive selection 

Diffractive events are characterised by a small momentum exchange at the proton vertex 
and by the presence of a large rapidity gap (LRG) between the proton beam direction and 
the hadronic hnal state. Diffractive DIS events were selected by the following additional 
criteria: 

• Xp < 0.01, where Xp is the fraction of the proton momentum carried by the diffract¬ 
ive exchange, calculated according to the formula xp = -|- Mx) / -|- W'^), in 

which Mx denotes the invariant mass of the hadronic state recoiling against the lead¬ 
ing proton and was reconstructed from the EFOs excluding the scattered electron 
candidate; this cut reduces the non-diffractive background; 

• hmax < 2, where rjras,^ is dehned as the pseudorapidity of the most forward EFO, with 
an energy greater than E-^yo = 400 MeV; this cut ensures the presence of a LRG in 
the event; 

• Mx > 5 GeV; this cut removes events with resonant particle production and ensures 
that there is enough energy in the system to create two jets with high transverse 
momenta. 

The origin of exclusive dijet events in diffraction is not unique. The most natural contribu¬ 
tion comes from exclusive production of quark-antiquark pairs, but other contributions, in 
particular from quark-antiquark-gluon, are not excluded. It is predicted im that the ratio 
of qq to qqg production changes signihcantly with the parameter 13 (or Mx) in contrast to 
other kinematic variables. To get insight into the origin of exclusive diffractive dijet events, 
the data were analysed as a function of /3, calculated according to /9 = {Q^ + Mx)- 
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4.3 Jet selection 


The /cT’-cluster algorithm known as the Durham jet algorithm [6Tl[62], as implemented in 
the Fast Jet package |63| . was used for jet reconstruction. Exclusive jets are of interest in 
this analysis, so the algorithm was used in “exclusive mode” i.e. each object representing 
a particle or a group of particles had to be finally associated to a jet. The algorithm is 
defined in the following way: first all objects were boosted to the 7 *-P rest frame. Then, 
the relative distance of each pair of objects, was calculated as 

krij = 2 min(£;|, £;|)(1 - cos %), 

where 9ij is the angle between objects i and j and Ei and Ej are the energies of the objects 
i and j. The minimum was found and if 

_ 

~ Ml 

objects i and j were merged. The merging of the 4-vectors was done using the recom¬ 
bination “E-scheme”, with simple 4-vector addition, which is the only Lorentz invariant 
scheme [M]. It causes the cluster objects to acquire mass and the total invariant mass, 
Mx, coincides with the invariant mass of the jet system. The clustering procedure was 
repeated until all yij values exceeded a given threshold, i/cut, and all the remaining objects 
were then labelled as jets. Applied in the centre-of-mass rest frame, this algorithm pro¬ 
duces at least two jets in every event. The same jet-search procedure was applied to the 
final-state hadrons for simulated events. 

Figure [3] shows the measured fractions for 2, 3 and 4 jets in the event as a function of 
the jet resolution parameter, ycut [H], in the region 0.01 < ycnt < 0.25. The rate of 
dijet reconstruction varies from 70% at ycut = 0.1 to 90% at ycut = 0.2. The measured 
jet fractions were compared to jet fractions predicted by SATRAP after reweighting of 
kinematic variables as described in Section [5l SATRAP provides a good description of the 
measurement. Jets were reconstructed with a resolution parameter fixed to ycut = 0.15. 
Events with exactly two reconstructed jets were selected. 

Finally, a lower limit of the jet transverse momenta in the centre-of-mass frame was re¬ 
quired, PTjet > 2 GeV. This value was chosen as a compromise between having a value of 
PTjet large enough so that perturbative calculations are still valid and on the other hand 
small enough so that a good statistical accuracy can be still obtained. 


5 Comparison between data and Monte Carlo 

Data and Monte Carlo predictions for several kinematic and jet variables were compared 
at the detector level. The MC event distributions which had been generated with SAT- 
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RAP were reweighted in a multidimensional space with respect to; inelasticity y, jet 
pseudorapidity ?7jet, PT,]et, (3 and x-p. In addition, the prediction of qq production 

from Bartels et ah j6] was used for reweighting in 0. 

The background originating from diffractive dijet photoproduction and non-diffractive 
dijet production was estimated from Monte Carlo simulations as described in Section [3l 
The background from beam-gas interactions and cosmic-ray events was investigated using 
data taken with empty proton-beam bunches and estimated to be negligible. 

For each of the distributions presented in this section, all selection criteria discussed above 
were applied except the cut on the shown quantity. The estimated background, normalised 
to the luminosity of this analysis, is also shown. 

Figure m shows the variables characterising the DIS events, y, W, E—Pz and .^vtx, 

while Fig. [5] shows the variables characterising the diffractive events, xp, Mx, (3 and prasa- 
Besides exclusive events, the data contains proton dissociation, e + p —)■ e-fjetl-|-jet2-|-y, 
for which the particles stemming from the process of dissociation disappear undetected in 
the proton beam hole. Except for r^max, the events with proton dissociation are expected 
to yield the same shape of the distributions as the exclusive dijet events, changing only 
the normalisation (according to the factorisation hypothesis, see Section [3]), and are not 
shown separately. These events were not treated as a background. The experimental dis¬ 
tributions were compared to the sum of the background distributions and the SATRAP 
MC. The background was normalised to the luminosity and the SATRAP MC to the num¬ 
ber of events remaining in the data after background subtraction. The r/max distribution 
(Fig. [5](d)) shows the distribution of events with proton dissociation which was determined 
separately as described in Section O In Fig. ISl(d), the r^max distribution is compared to the 
normalised sum of three contributions including that of events with proton dissociation. 
All data distributions, except for y and Pra^xi are reasonably well described by the MC 
predictions. Most of the difference between data and MC in the y distribution is outside 
the analysed region {y > 0.64). The incorrect y description at the cut value was taken 
into account in the systematic uncertainty which was determined by varying the cut. The 
shift in the 77 m ax distribution is accounted for in the systematic uncertainty of the proton 
dissociation background (see Section [6]). 

In Fig. ini jet properties in the 7 *-P centre-of-mass system are presented: the distributions 
of the jet angles 9 and 0, the number of EFOs clustered into the jets and the jet transverse 
momentum PTjet- All distributions are reasonably well described by the sum of SATRAP 
events and the background distribution. The difference between data and MC for values 
of 0 close to 0 is not expected to affect the result of unfolding in this quantity (see 
Section [7]). 

Jets reconstructed in the 7 *-P rest frame were transformed back to the laboratory (LAB) 
system. In Fig. [7|the distributions of the jet pseudorapidity and the jet transverse energy 
are shown in the laboratory system separately for higher- and lower-energy jets. They are 
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well described by the predicted shape. 


The jet algorithm used allows the association of the individual hadrons with a unique 
jet on an event-by-event basis. To study the topology of the jets, the energy flow of 
particles around the jet axes was considered in both the centre-of-mass system and the 
laboratory system. In this study, Arj and Aip denote the differences between the jet 
axis and, respectively, the pseudorapidity and the azimuthal angle of the EFOs in the 
event. In Figs. [8] the energy flows around the axis of the reference jet, that is the jet with 
positive Z-component of the momentum, are shown in the 7 *-P centre-of-mass system. 
The corresponding distributions in the laboratory system are presented in Figs. [9l It 
is observed that energy flows around the reference jet axis are well reproduced by the 
SATRAP MC. As expected, the jets are produced back-to-back in the 7 *-lP centre-of- 
mass system, and are quite broad. However, in the laboratory system, most of their 
energy is concentrated within a cone of radius approximately equal to one unit in the rj-ip 
plane with distance dehned as r = Arf -|- A(p^. 

The quality of the description of the data by the MC gives conhdence in the use of the 
MC for unfolding differential cross sections to the hadron level (see Section [7]). 


6 Estimate of dijet production with proton dissociation 

The contribution of events with a detected dissociated proton system is highly suppressed 
due to the nominal selection cuts applied to the data, i.e. by requiring exactly two jets, 
X]p < 0.01 and ?7max < 2, and has been considered to be negligible. However, the con¬ 
tribution of proton-dissociative events, where the proton-dissociative system escapes un¬ 
detected, is not negligible. It was estimated using EPSOFT, after further tuning of the 
distribution of the mass of the dissociated proton system. My. The simulation shows 
that due to the acceptance of the calorimeter, determined by the detector geometry, a 
dissociated proton system of mass smaller than about 6 GeV stays undetected. 

In order to estimate the amount of dissociated proton events, a sample enriched in such 
events was selected as follows. Kinematic variables and jets were reconstructed from the 
EFOs in the range rj < 2. All selection cuts described in Section 01 except the r^max cut, 
were then applied. In order to suppress non-diffractive contributions to the dijet sample, 
events with EFOs in the range 2 < rj < 3.5 were rejected. The remaining sample of 
events with EFOs in the range rj > 3.5 consisted almost entirely of diffractive dijets with a 
detected dissociated proton system. From the comparison of the energy sum of all EFOs 
with T] > 3.5 between data and simulated events, the following parameterisation of the 
My distribution was extracted: 

d(J'yp_^jetijet2 V f 

dM| Mp' 





The fraction of simulated events with proton dissociation was determined by a fit to the 
distribution of r/max shown in Fig. O^d). 

The systematic uncertainty of this fraction was estimated in the following steps: 

• the shape of the My distribution was varied by changing the exponent by ±0.6, 
because in this way the of fhe comparison between data and EPSOFT simulation 
was raised by 1; 

• the fit of the fraction was repeated taking into account a shift of r^max by ±0.1 
according to the observed shift between data and simulated events. 

Both uncertainties were added in quadrature. 

The fraction of events with r/max < 2 associated to the proton-dissociative system, which es¬ 
caped undetected in the beam hole, was estimated to be /pdiss = 45%±4%(stat.)±15%(syst.) 
No evidence was found that /pdiss depends on or jS. Therefore, in the following sections, 
the selected data sample was scaled by a constant factor correcting for proton-dissociative 
events. 


7 Unfolding of the hadron-level cross section 

An unfolding method was used to obtain hadron-level differential cross sections for pro¬ 
duction of dijets, reconstructed with jet-resolution parameter i/cnt = 0.15, as a function of 
jS and (j) in the following kinematic region: 

• Q2 > 25 GeV^; 

• 90 < IT < 250 GeV; 

• Xp < 0.01; 

• Mx > 5 GeV; 

• Ajets = 2; 

• Prjet > 2 GeV. 

The unfolding was performed by calculating a detector response matrix, which represents 
a linear transformation of the hadron-level two-dimensional distribution of or f3- 

PTjet to a detector-level distribution. The response matrix was based on the weighted 
SATRAP MG simulation. It includes effects of limited detector and trigger efficiencies, 
hnite detector resolutions, migrations from outside the phase space and distortions due 
to QED radiation. The unfolding procedure was based on the regularised inversion of the 
response matrix using Singular Value Decomposition (SVD) as implemented in the TS- 
VDUnfold package [65]. The implementation was prepared for one-dimensional problems 
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and the studied two-dimensional distributions were transformed into one-dimensional dis¬ 
tributions [nn]- The regularisation parameter was determined according to the procedure 
suggested by the authors of the unfolding package. 

The used unfolding method takes into account the imperfect hadron level MC simulation 
and corrects for it. 


8 Systematic uncertainties 


The systematic uncertainties of the cross sections were estimated by calculating the dif¬ 
ference between results obtained with standard and varied settings for each bin of the 
unfolded distribution, except for the uncertainty on /pdiss, which was assumed to give a 
common normalisation uncertainty in all the bins. 

The sources of systematic uncertainty were divided into two types. Those originating 
from detector simulation were investigated by introducing changes only to MC samples 
at the detector level, while the data samples were not altered. The following checks were 
performed: 

• the energy scale of the calorimeter objects associated with the jet with the highest 
transverse momentum in the laboratory frame was varied by ±5%; the corresponding 
systematic uncertainty is in the range of +2%, —8%; 

• the jet transverse momentum resolution was varied by ±1%, because in this way 
the of fho comparison of data to MC in the distribution of the jet transverse 
momentum was raised by 1. 

Systematic effects originating from event-select ion cuts were investigated by varying the 
criteria used to select events for both data and simulated events in the following ways: 

• Q2 > 25 ± 1.7 GeV^; 

• 90 ± 7.4 GeV < kk < 250 ± 8.4 GeV; 

• 0.1 ±0.04 <y< 0.64 ±0.03; 

• I 2’vtx I < 30 ± 5 cm; 

• Xp < 0.01 ± 0.001; 

• ^max ^ 2 ± 0.2, 

• Mx > 5 ± 0.8 GeV; 

• > 0.4 ±0.1 GeV. 


10 



The uncertainty related to the M^-cut variation is in the range of ±5%. The total un¬ 
certainty related to the event-selection cuts excluding the Mx cut is smaller than ±6%. 
Uncertainties originating from the most significant sources are presented in Fig. [101 

Positive and negative uncertainties were separately added in quadrature. The correspond¬ 
ing total systematic uncertainty is also shown in Fig. [101 The normalisation uncertainty of 
the cross section related to the luminosity (see Section [2]) as well as to /pdiss is not shown 
on the following figures but is included as a separate column in the tables of cross sections. 
The total uncertainties of the measured cross sections are dominated by the systematic 
component. 


9 Cross sections 

Cross sections were measured at the hadron level in the kinematic range described in Sec¬ 
tion [71 Backgrounds from diffractive photoproduction and non-diffractive dijet production 
were subtracted. 

In order to calculate the cross sections for exclusive dijet production, the measured cross 
sections were scaled by a factor of (1 — /pdiss) = 0.55 according to the estimate of the 
proton-dissociative background described in Section [6l 

The values of the cross-sections dcr/d/3 and da/d/) in five bins of jS are given in Tables [T] 
and [21 and shown in Fig. [TTl The statistical uncertainties presented in the figures cor¬ 
respond to the diagonal elements of the covariance matrices, which are available in elec¬ 
tronic format [66]. The dcr/d/3 distribution is, due to the kinematics, restricted to the 
range 0.04 < (3 < 0.92. The da/dcf) distribution is shown in five bins of f3 in the range 
0.04 < (3 < 0.7. The cut at 0.7 excludes a region with a low number of events. 

The (j) distributions show a significant feature: when going from small to large values of 
/9, the shape varies and the slope of the angular distribution changes sign. The variation 
of the shape was quantified by fitting a function to the 0 distributions including the 
full statistical covariance matrix and the systematic uncertainties, the latter by using 
the profile method EZj. The fitted function is predicted by theoretical calculations (see 
Section H]) to be proportional to (1 -|- A cos 20). The data are well described by the fitted 
function. The resulting values of A are shown in Tableland Fig. [T2J The parameter A 
decreases with increasing (3 and changes sign around (3 = 0.4. 

10 Comparison with model predictions 


The differential cross sections were compared to MC predictions for the Resolved-Pomeron 
model and the Two-Gluon-Exchange model. In the Resolved-Pomeron model [18], the 
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diffractive scattering is factorised into a Ponieron flux from the proton and the hard 
interaction between the virtual photon and a constituent parton of the Pomeron. An 
example of such a process is shown in Fig. [131 where a qq pair is produced by a boson- 
gluon fusion (BGF) process associated with the emission of a Pomeron remnant. This 
model requires the proton diffractive gluon density as an input for the calculation of the 
cross section. The predictions considered in this article are based on the parameterisation 
of the diffractive gluon density obtained from hts (HI 2006 hts A and B) to HI inclusive 
diffractive data [13]. The shape of the 0 distribution is essentially identical in all models 
based on the BGF process, including both the Resolved-Pomeron and the Soft Golour 
Interactions (SGI) model |68| . 

In the Two-Gluon-Exchange model |HI61IT7]. the diffractive production of a qq pair is due 
to the exchange of a two-gluon colour-singlet state. The process is schematically shown in 
Fig. [m The qq pair hadronises into a dijet hnal state. For large diffractive masses, i.e. 
at low values of /3, the cross section for the production of a qq pair with an extra gluon 
is larger than that of the qq production. The diagram of this process is shown in Fig. [151 
The qqg hnal state also contributes to the dijet event sample if two of the partons are not 
resolved by the jet algorithm (see Section KM- 


The qq pair production was calculated to second order in QGD, using the running strong- 
interaction coupling constant ^^(/i) with the scale p = [HE], where pt 

denotes the transverse momentum of the quarks in the 7 *-IP rest frame with respect to 
the virtual photon momentum and Mx is the invariant mass of the diffractive system. 
The prediction was calculated with a cut on the transverse momentum of the quarks, 
Pt > 1 GeV. The cross section is proportional to the square of the gluon density of the 
proton, g{x-p,pi^) 


da oc 


Pt 


The cross section for the qqg hnal state was calculated taking into account that it is 
proportional to the square of the gluon density of the proton, p(a;p,/c|n), at a scale 
which effectively involves the transverse momenta of all three partons. For the calculation 
of the cross section, a hxed value of ag = 0.25 El and the GRV [61] parameterisation of 
the gluon density were used and the same cut was applied on the transverse momentum 
of all partons: Pt, parton > Pt, cut with the value adjusted to the data (see Section flU.ip . 
In contrast to qq production, the exclusive dijet cross section calculated for the qqg hnal 
state is sensitive to the parton-level cut pT,cut- This is a consequence of the fact that two 
of the partons form a single jet. 
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10.1 Contribution of the qq dijet component in the prediction of 
the Two-Gluon-Exchange model 

In the Two-Gluon-Exchange model, the 0 distribution predicted for qq and qqg have dif¬ 
ferent shapes. This allows the ratio Rqq = o'{qq)/a{qq + qqg) to be determined by studying 
the measured 0 distributions. The results are shown in Fig. [TSl The ratio was measured 
only in the region of /3 G (0.3, 0.7) since elsewhere the uncertainty estimation is unreliable 
due to the measured value being too close to 0 or 1. The ratio Rqq predicted by the model 
depends on the parton transverse-momentum cut applied. The Pt, cut value of \/2 GeV 
used in the original calculation mi significantly underestimates the ratio. A scan of the 
parton transverse-momentum cut showed that the measured ratio can be well described 
throughout the considered range with pT,cut = 1-75 GeV. Both this value of pt, cut and the 
original value were used for calculating the Two-Gluon-Exchange model predictions. 

10.2 Differential cross-section dcr/d/3 

The cross-section dcr/d/3 is shown in Fig. [T7] together with the predictions from both 
models. The prediction of the Resolved-Pomeron model decreases with increasing f3 faster 
than the measured cross section, for both fit A and fit B. The difference between data 
and prediction is less pronounced for fit A than for fit B, which is consistent with the 
observation that the ratio of gluon densities increases with increasing [3 |13]. Predictions 
and data differ by a factor of two for small values of f3 and about ten for large values. 

The Two-Gluon-Exchange model prediction, which includes qq and qqg, describes the 
shape of the measured (3 distribution reasonably well. The predicted integrated cross sec¬ 
tion is (T = 38 pb, while the measured cross section is a = 72 pb with a normalisation uncer¬ 
tainty originating from the proton-dissociation background of M(/pdiss)/(l — /pdiss) = 27%, 
where ^(/pdiss) is fhe uncertainty in the fraction of events with a dissociated proton. Al¬ 
though the difference between the predicted and measured cross section is not significant, 
it could indicate that the NLO corrections are large or the cross-section enhancement 
arising from the evolution of the off-diagonal gluon distribution is significant [7]. The 
prediction based on qq production alone fails to describe the shape of the distribution at 
low values of (3 but is almost sufficient to describe it at large (3, where the qqg component 
is less important. 

10.3 Differential cross-section dcr/d(f) 

The cross-sections da/dcf) are shown in Fig. [T7] in five different (3 ranges together with 
the predictions of both models. The comparison of the shapes has been quantified by 
calculating the slope parameter A. The results are shown in Fig. [THl The Resolved- 
Ponieron model predicts an almost constant, positive value of A in the whole f3 range. The 
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Two-Gluon-Exchange model {qq -|- qqg) predicts a value of A which varies from positive 
to negative. In contrast to the Resolved-Pomeron model, the Two-Gluon-Exchange model 
agrees quantitatively with the data in the range 0.3 < f3 < 0.7. The prediction based 
on qq production alone describes the shape of the distributions at large /9, where the qqg 
component is less important. 


11 Summary 

The hrst measurement of diffractive production of exclusive dijets in deep inelastic scatter¬ 
ing, 'j*+p —jetl-|-jet2-|-p, was presented. The differential cross-sections da/d/S and dcr/d0 
in bins of f3 were measured in the kinematic range: > 25 GeV^, 90 < IT < 250 GeV, 

Mx > 5 GeV, xp < 0.01 and prjet > 2 GeV using an integrated luminosity of 372 pb“^. 

The measured absolute cross sections are larger than those predicted by both the Resolved- 
Pomeron and the Two-Gluon-Exchange models. The difference between the data and 
the Resolved-Pomeron model sX (3 > 0.4 is significant. The Two-Gluon-Exchange model 
predictions agree with the data within the experimental uncertainty and are themselves 
subject to possible large theoretical uncertainties. The shape of the 0 distributions was 
parameterised with the function 1 + A cos 20, as motivated by theory. The Two-Gluon- 
Exchange model predicts reasonably well the measured value of A as a function of 0, 
whereas the Resolved-Pomeron model exhibits a different trend. 
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Table 1: Differential cross-section da/dfi in the kinematic range: > 25GeV^, 

90 < IT < 250 GeV, xp < 0.01, Mx > 5 GeV and Prjet > 2 GeV. The statistical uncer¬ 
tainties are given by the diagonal part of the covariance matrix. Systematic uncertainties 
are explained in the text. The contribution from proton dissociation was subtracted. The 
uncertainty of the subtraction determines the uncertainty of the normalisation also given 
in the table. 


f3 da/d/S (pb) 


0.04 - 

0.15 

159.7±1.8(stat.) 

l4;7(sys.)±45.2(norm.) 

0.15 - 

0.3 

175.1±1.3(stat.) 

lg;o(sys.)±49.6(norm.) 

0.3 - 

0.4 

132.3±1.2(stat.) 

lg;5’(sys.)±37.5(norm.) 

0.4 - 

0.5 

82.1±0.9(stat.) 

l5;o(sys.)±23.3(norm.) 

0.5 - 

0.7 

29.0±0.5(stat.) 

l|2(sys-)=t8-3(norm.) 

0.7 - 

0.92 

2.47±0.06(stat.; 

1 l[5;^?(sys.)±0.70(norm.) 
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Table 2: Differential cross-section da/df in the kinematic range: > 25GeV^, 

90 < IT < 250 GeV, xp < 0.01, Mx > 5 GeV and pT,jet > 2 GeV. Statistical uncertain¬ 
ties are given by the diagonal elements of the covariance matrix. Systematic uncertainties 
are explained in the text. The contribution from proton dissociation was subtracted. The 
uncertainty of the subtraction determines the uncertainty of the normalisation given in the 
table. 


0 

(rad) 

dcr/d 0 (pb/rad) 

0.04 <13 < 0.15 

0 

0.314 

0.628 

0.942 

1.26 

- 0.314 

- 0.628 

- 0.942 

- 1.26 

- 1.57 

14.64±0.64(stat.) ^Q;5Q(sys.)±4.15(norm.) 
12.73±0.49(stat.) lQ;gi(sys.)±3.60(norm.) 
10.71±0.43(stat.) ^o;82(sys.)±3.03(norm.) 
9.46±0.39(stat.) tQ; 53 (sys.)± 2 . 68 (norm.) 
8.89±0.45(stat.) lQ;45(sys.)±2.52(norm.) 

0.15 < /3 < 0.3 

0 

0.314 

0.628 

0.942 

1.26 

- 0.314 

- 0.628 

- 0.942 

- 1.26 

- 1.57 

21.03±0.60(stat.) t|; 43 (sys.)± 5 . 95 (norm.) 
17.01±0.44(stat.) ^j;4g(sys.)±4.82(norm.) 
14.89±0.41(stat.)l|);|j[((sys.)±4.22(norm.) 
15.20±0.39(stat.) ^o;7g(sys.)±4.30(norm.) 
15.33±0.49(stat.) to;84(sys.)±4.34(norm.) 

0.3 < /3 < 0.4 

0 

0.314 

0.628 

0.942 

1.26 

- 0.314 

- 0.628 

- 0.942 

- 1.26 

- 1.57 

9.61±0.43(stat.) lo;84(sys.)±2.72(norm.) 
8.18±0.29(stat.) lQ;7Q(sys.)±2.32(norm.) 
7.78±0.28(stat.) ^o;58(sys.)±2.20(norm.) 
8.36±0.29(stat.) Io;g 3 (sys.)± 2 . 37 (norm.) 
8.39±0.41(stat.) to;77(sys.)±2.38(norm.) 

0.4 < /3 < 0.5 

0 

0.314 

0.628 

0.942 

1.26 

- 0.314 

- 0.628 

- 0.942 

- 1.26 

- 1.57 

5.68±0.33(stat.) ^Q;g 7 (sys.)±1.61 (norm.) 
4.60±0.23(stat.) tQ;g4(sys.)±1.31(norm.) 
4.67±0.22(stat.) lQ; 5 g(sys.)± 1.33(norm.) 
5.64±0.24(stat.) ^Q;g 7 (sys.)± 1.60(norm.) 
5.32±0.36(stat.) lQ;gQ(sys.)±1.51(norm.) 

0.5 < /3 < 0.7 

0 

0.314 

0.628 

0.942 

1.26 

- 0.314 

- 0.628 

- 0.942 

- 1.26 

- 1.57 

3.58±0.22(stat.) lQ;gg(sys.)±1.02(norm.) 
2.81±0.16(stat.) ^Q; 55 (sys.)± 0 . 80 (norni.) 
3.20±0.17(stat.) ^Q;54(sys.)±0.91(norm.) 
4.19±0.19(stat.) lQ; 54 (sys.)± 1.19(norm.) 
4.78±0.28(stat.) ^Q; 53 (sys.)± 1.36(norm.) 


20 















Table 3: Results of the fit to the cross-section da/dfi in bins of fi. The fitted function is 
proportional to {1 + A cos 20). The uncertainty includes both statistical and systematical 
contributions (see text). 




0.04 - 

0.15 

0.15 - 

0.3 

0.3 - 

0.4 

0.4 - 

0.5 

0.5 - 

0.7 


A 

0.256±0.030 

0.130±0.028 

0.053±0.045 

-0.037±0.054 

-0.196±0.070 
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Figure 1: Schematic view of the diffractive production of exclusive dijets in electron- 
proton DIS. 



Figure 2: Definition of planes and angles in the 7 *-P centre-of-mass system. The lepton 
plane is defined by the 7 * and e momenta. The jet plane is defined by the 7 * and dijet 
directions. The angle (j) is the angle between these two planes. The jet polar angle, 9, is 
the angle between the directions of the jets and 7 *. 
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Figure 3: The probability of finding two, three and four jets in the final state as a function 
of the ycut parameter (see text). Data are shown as full dots. Statistical errors are smaller 
than the dot size. Predictions of SATRAP are shown as histograms. The distributions are 
not corrected for detector effects. 
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Figure 4: Comparison between data (dots) and the sum (line) of the SATRAP MC and 
background contributions (shaded), where events with a dissociated proton are not treated 
as background, for kinematic variables: (a) exchanged photon virtuality, , (b) scattered 
electron energy, E(, (c) inelasticity, y, (d) invariant mass of the 7 *^ system, W, (e) 
the guantity E — Pz and (f) the Z-coordinate of the interaction vertex. The error bars 
represent statistical errors (generally not visible). The background was normalised to the 
luminosity and the SATRAP MC to the number of events remaining in the data after 
background subtraction. All selection cuts are applied except for the cut on the variable 
shown in each plot. 


24 

































ZEUS 





w 8000 



11 

'max 


Figure 5: Comparison between data (dots) and the sum (line) of the SATRAP MC and 
background contributions (dark histogram), where events with a dissociated proton were 
not treated as background, for the kinematic variables of the diffractive process: (a) frac¬ 
tion loss in the longitudinal momentum of the proton, xp, (b) the invariant mass of the 
diffractive system, Mx, (c) the Bjorken-like variable, (3, and (d) the pseudorapidity of the 
most forward EFO, In (d), the component from events with a dissociated proton is 

also shown (light histogram). Other details as for Fig. 
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Figure 6: Comparison between data and MC for variables characterising the jet properties 
in the 7 *-P rest frame: (a) the azimuthal angle, (f, (b) the polar angle, 9, (c) the number 
of EFOs clustered into a jet and (d) the jet transverse momentum, PTjet- Other details as 
for Fig. g} 
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Figure 7: Comparison between data and MC for variables characterising the properties of 
the higher and the lower energy jets in the LAB frame: (a) and (b) the jet pseudorapidity, 
(b) and (c) the jet transverse energy. Other details as for Fig. 
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Figure 8: The energy flow in the 7 *-P rest frame around the jet axis, averaged over all 
selected dijet events, is shown as a function of distances in azimuthal angle and pseu¬ 
dorapidity (Aip and Ap). In both cases, the energy flow is integrated over the full available 
range of the other variable. Data for both jets are shown as full dots. Statistical uncertain¬ 
ties are smaller than point markers. The energy flow of EFOs belonging to the reference jet 
only are shown as full sguares, where the reference jet was chosen as the jet with positive 
Pz momentum. Predictions of the SATRAP MC are shown as histograms. 
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Figure 9: The energy flow in the laboratory frame, around the jet axis, averaged over 
all selected dijet events, is shown as a function of distances in azimuthal angle and pseu¬ 
dorapidity (Aip and Arj). In both cases, the energy flow is integrated over the full available 
range of the other variable. Data for both jets are shown as full dots. Statistical uncertain¬ 
ties are smaller than point markers. The energy flow of EFOs belonging to the reference jet 
only are shown as full squares, where the reference jet was chosen as the jet with positive 
Pz momentum. Predictions of the SATRAP MC are shown as histograms. 
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Figure 10: The most significant sources of systematic uncertainties for da/dfi and da/df 
in five fd ranges. Total systematic uncertainties and total uncertainties (systematic and 
statistical added in quadrature) are shown as shaded and dark-shaded bands. 
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Figure 11: Differential cross sections for exclusive dijet production: da/ dfd (in log scale) 
and da/df (in linear scale) in five bins of fi. Contributions from proton-dissociative dijet 
production were subtracted. The full line represents the fitted function proportional to 
1 + A cos 20. Statistical and systematic uncertainties were included in the fit. The total 
error bars show statistical and systematic uncertainties added in quadrature. The statist¬ 
ical uncertainties were taken from the diagonal elements of the covariance matrix. The 
systematic uncertainties do not include the uncertainty of the subtraction of the proton- 
dissociative contribution. This normalisation uncertainty is shown as a grey band only in 
the da/dfi distribution. 
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Figure 12: The shape parameter A as a function of 13 resulting from the fits to da/dcj) 
with a function proportional to 1 + A cos 20. The statistical and systematic uncertainties 
were included in the fit. 




Figure 13: Diagram of diffractive boson-gluon fusion in the Resolved-Pomeron model. 
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Figure 14: Example diagram of qq production in the Two-Gluon-Exchange model. 
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Figure 16: The Rqg = (^{qq)/{a{qq) + (y{qqg)), determined in a fit of the predicted shapes 
to the measured distributions given in Fig. [771 The fit takes into account the full covari¬ 
ance matrix. The predicted ratio is shown for two choices of pT,cut-' for the \/2 GeV used 
for the published calculations mi and for 1.75 GeV, determined in a fit. 
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Figure 17: Differential cross sections as in Fig. [77] in comparison to model predictions 
da/df (in log scale) and da/dcj) in bins of (3 (in linear scale). Contributions from proton- 
dissociative dijet production were subtracted. The systematic uncertainties do not include 
the uncertainty due to the subtraction. The Two-Gluon-Exchange model is presented with 
PT,cut = T75 GeV. The bands on theoretical expectations represent statistical uncertainties 
only. 


35 








































































ZEUS 


< 0.4 
0.2 
0 

- 0.2 

- 0.4 

Figure 18: The shape parameter A as a function of f3 in comparison to the values of A 
obtained from distributions predicted by the Resolved-Pomeron model and the Two-Gluon- 
Exchange model. The bands on BGF Fit B and two-gluon pT,cut = T75 GeV represent 
statistical uncertainties. 
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